Effective short-range interaction for spin-singlet P-wave nucleon-nucleon scattering 
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Distorted-wave methods are used to remove the effects of one- and two-pion exchange up to order 
Q 3 from the empirical x Pi phase shift. The one divergence that arises can be renormalised using 
an order-Q 2 counterterm which is provided by the (Weinberg) power counting appropriate to the 
effective field theory for this channel. The residual interaction is used to estimate the scale of the 
underlying physics. 



Effective field theories (EFTs) provide important tools for constructing nuclear forces within a systematic frame- 
work. 1 They respect the symmetries of underlying theory, in this case Quantum Chromodynamics, and offer model- 
independent descriptions of the dynamics at low energy scales. They rely on the existence of clear separation of scales 
between the low-energy physics of interest and the underlying, high-energy (or short-distance) physics. This makes it 
possible to use a perturbation theory where the expansion parameters are ratios of low- to high-energy scales. 

To construct such a theory one begins by writing down the most general Lagrangian or Hamiltonian for the relevant 
low-energy degrees of freedom, containing all the interaction terms allowed by the symmetries of the system. The 
parameters of the theory are the (initially unknown) constants multiplying these terms. There are in general an 
infinite number of these but they can be organised in terms of a "power counting" of low-energy scales. To any given 
order in this counting, only a finite number of terms are needed. 

Since Weinberg first suggested that these ideas could be applied to nuclear forces @], there has been considerable 
debate about the appropriate power counting to use. 2 Most of these arguments have been about channels where 
scattering is strong and some pieces of the potential may need to be treated nonperturbatively, namely the S waves 
and the lower spin-triplet waves. In the higher partial waves, especially the spin-singlet ones, the scattering is weak 
and the picture is clearer. For these channels, Weinberg's original power counting is valid. This is based on simple 
dimensional analysis, counting powers of momenta and the pion mass. These low-energy scales are generically denoted 
byQ. 

Many applications of EFTs test convergence by presenting theoretical nucleon-nucleon phase shifts, calculated at 
leading and higher orders, alongside empirical data. At each order, short-range interactions, expanded as polynomials 
in energy or momenta, have to be refitted to the data. A powerful alternative is to use distorted-wave (DW) methods to 
remove the effects of known long-range forces from the observed phase shifts to leave residual short-range interactions. 
The sizes and energy-dependences of these can then be used to test the convergence of the EFT and to estimate the 
scales of the underlying physics. Such an approach was proposed in Ref. [7J and applied to the peripheral spin-singlet 
waves. More recently it has been applied to spin-triplet waves [8] and the 1 5o channel [9] and closely related methods 
have been used by other groups to study the S waves [TU1 E] ■ 

In this note, we extend the method of Ref. [7J by applying it to the 1 P\ channel. We first construct DWs for the 
leading-order one-pion-exchange (OPE) potential. Iterating this potential to all orders is not essential in a spin-singlet 
wave with L ^ since Weinberg's power counting is expected to hold. However doing so avoids the need to calculate 
terms to fourth order in perturbation theory. Since only the central part of OPE contributes and its 1/r singularity 
does not alter the power-law behaviour of the wave functions near the origin, this treatment does not affect the power 
counting for the short-range pieces of the potential. 

From the K matrix that describes scattering between these DWs, we can define a residual interaction. This still 
contains effects arising from long-range potentials of orders Q 2 and above. We then use the distorted-wave Born 
approximation (DWBA) to subtract the matrix elements of the order-Q 2 and Q 3 two-pion-exchange (TPE) potentials 
[T^HH]. This leaves a residual interaction that should represent purely short-range physics to this order. 

One issue that arises here, but not for the channels studied in Ref. [7j, is that the matrix elements of TPE between 
the DWs are divergent. This divergence must be removed by renormalisation so that only finite quantities are being 
treated in perturbation theory. Weinberg's power counting provides one energy-independent counterterm, of order 
Q 2 , in this channel and this is sufficient to cancel the divergence. This is similar to what has been found in other 



1 For reviews, see Refs. 

2 For recent summaries of two of these points of view, see Refs. [5] and [6]. 
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waves where different counting schemes apply [8HTT]. 3 

We start by outlining the main features of the method from Ref. [7]. The DWs ip OPE (p, r) are obtained by solving 
the radial Schrodinger equation with the OPE and centrifugal potentials, 

d V'ope 2 d-ip 

OPE 

/ L(L + 1) (o) 



dr 2 r dr 



( L(L r2 +1) + M N V^ E {r^j ^ OPE (p,r) = p 2 ^ OPE (p,r). (1) 



Here L = 1 is the orbital angular momentum, M N is the nucleon mass and p is the on-shell relative momentum in the 
centre-of-mass frame. We write the leading-order OPE potential in the same form as used in the Nijmegen analyses 
and take their preferred value for the irN coupling, f 2 NN = 0.075 [T71H5] . 

From the large-r forms of these waves we can extract the OPE phase shift, 5 OPE (p). Taking the difference between 
this and the empirical phase shift, S(p), we can define a residual iC-matrix, 

An 

K{ ?) = ~TT^ tan W - *op*(p)). (2) 

lvl N p 

This describes the additional scattering between the DWs, produced by short-range interactions and long-range forces 
of order Q 2 and higher. If we take a (5-shell form for the potential responsible 

V s (p,r) = ^V(p)S(r-R ), (3) 
then we can extract its strength directly from K(p), 



V {2) (P) = QI J" R , 2 K(P), (4) 



Ro 

,R ) 2 

where the superscript (2) indicates that long-range effects of order Q 2 and higher are still present. 

Finally, we use the DWBA to subtract the effects of long-range potentials to order Q 3 . This leaves a residual 
short-range potential whose strength is given by 

f(4) (P) = 51 f p ^ ( R (P) ~ (^op-CpJIV^ + Vt ( p E 3 W e (p)>) , (5) 
9yor E (P, Ro) V ' 

where the order-Q 2 ' 3 TPE potentials can be found in Refs. jT2l [13] and the corresponding order-Q 2 recoil correction 
to OPE is given by Friar [T3] . 

The TPE potential has terms with 1/r 6 and 1/r 5 singularities. Consequently, the radial integrals in its matrix 
elements between P waves, which have the form ip(r) oc r for small r, contain 1/r and lnr divergences. To rectify 
this, we regularise the integrals by imposing a radial cut-off at the same value of Rq as used in the <5-shell potential. 
We can then renormalise the matrix elements by subtracting the sum of these divergent pieces (a constant which is 
independent of energy). 
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FIG. 1: The short-range effective potential V^ 2 \p), in fm 4 , plotted against the lab kinetic energy T in MeV. The four curves 
correspond to different Nijmegen PWAs. A cut-off radius of -Ro = 0.1 fm was used. 



3 Note that here we follow the power counting by treating the order-Q 2,3 potentials as perturbations. Other approaches, which treat the 
whole potential to all orders, can lead to different conclusions 15. [161 . 
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We show first, in Fig. [I] the effective short-range potential V^(p) in the 1 P\ partial wave after the removal of 
OPE only. The cut-off radius of the 5-shell was taken to be Rq — 0.1 fm. Results are shown for PWA93 and three 
different Nijmegen potentials [T7HT5] . all of which provide high-quality fits to the NN scattering data. One can see 
that significant energy dependence is present at low energies, like that found in the higher partial waves [Tj- This 
suggests that there may still be important contributions from long-range forces in this interaction. However we 
should note that the different Nijmegen analyses diverge for energies below about 50 MeV and so this region should 
not be regarded as well constrained by the data. This is because the centrifugal barrier prevents the nucleons from 
approaching each other closely at these energies and so the 1 P\ phase shift is completely dominated by OPE. 




FIG. 2: The renormalised short-range effective potential V^\p), in fm 4 , plotted against the lab kinetic energy T in MeV. 
The different curves show the potentials obtained with cut-off radii Ro = 0.8 fm (dotted), 0.4 fm, 0.2 fm, 0.1 fm, and 0.05 fm 
(solid). 

Using the DWBA to remove the effects of the order-Q 2,3 long-range potentials, we obtain the results shown in 
Figs. [2] and [3] The matrix elements of the TPE potential have been renormalised by simply subtracting their values 
at some very low energy, T = 5 MeV. Fig. [2] shows that the renormalised potential does indeed converge to a cut-off- 
independent form as Rq —> 0. For numerical convenience we present our results for R = 0.1 fm, for which any cut-off 
artefacts (terms in the potential proportional to positive powers of Ro) are very small. 




FIG. 3: The short-range effective potential V fi \p) for four different Nijmegen PWAs. A cut-off radius of Ro = 0.1 fm was 
used. 

The resulting renormalised potentials (j>) obtained from several Nijmegen PWAs are shown in Fig. [3] The 
"turnover" at low energies that was visible in Fig. [T] has been significantly reduced by subtracting the higher-order 
long-range forces. However we should point out that much of this effect lies in the region where uncertainties in the 
PWAs make it hard to draw very strong conclusions. 

To try to quantify the energy dependence of the residual potential, we made a least squares fit of a quadratic in 
the energy to V^(p). We fitted this over the energy range T — 100 — 200 MeV, which was chosen to avoid both 
the low-energy region where the data are too inaccurate to determine the short-range potential reliably, and the 
high-energy region where our EFT is expected to converge slowly if at all. A quadratic was chosen because attempts 
to fit higher order polynomials to this range did not lead to stable values for the coefficients. 

Writing the polynomial in the form 



V^(p) =a + a lP 2 + a 2 p 4 , 



(6) 
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we obtain ao = 0.28 fm 4 , a-y = —0.20 fm 6 and a-i — 0.0056 fm 8 . This is dominated by a linear dependence on the 
energy (p 2 ), consistent with the visual impression of Fig. [3] Without knowing the uncertainties on the phase shifts 
that are used to determine the potential, it is difficult to assign definite errors to these coefficients. However simple 
estimates suggest that the values for the first two are reasonably accurate but 02 is not. 

For present purposes, our main interest in these coefficients is to estimate the scale of the short-distance physics 
represented by V^(p). This is the scale that controls the convergence of the expansion of our EFT. Here, we cannot 
use ao for this because it depends on our (somewhat ad hoc) choice of renormalisation scheme. Also, as just noted, 
the coefficient 02 is not accurately determined and so we are left with only a±, the coefficient of the term linear in 
energy. 

If our theory has a "natural" scale dependence, ai should be of the form a\ — ai/Ag, where Ao is the scale of the 
underlying physics and Si is a dimensionless number of order unity. The higher-order pion-exchange forces that have 
not been subtracted out can contribute to this coefficient but, since these are at least of order Q , their contributions 
are suppressed by (m^/Ao) 2 . If we set a\ = —I, then our value of a\ implies a scale of approximately 260 MeV. This 
is similar to the breakdown scale that was recently estimated for the 1 So wave [9]. Such low values for the scales 
suggest the presence of additional low-energy physics that has not been considered in the present EFT. One important 
possible example of this is the A-resonance. An extended theory, including the A as an explicit degree of freedom as 
in Ref. [20], may lead to more natural short-range forces. 
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